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Abstract 



m A version of the virial theorem, which takes into account the effects of the non-compact extra- 

■ dimensions, is derived in the framework of the brane world models. In the braneworld scenario, 

the four dimensional effective Einstein equation has some extra terms, called dark radiation and 
dark pressure, respectively, which arise from the embedding of the 3-brane in the bulk. To derive 
the generalized virial theorem we use a method based on the collisionless Boltzmann equation. 



, The dark radiation term generates an equivalent mass term (the dark mass), which gives an ef- 

^Sj , fective contribution to the gravitational energy. This term may account for the well-known virial 
> 

00 \ theorem mass discrepancy in actual clusters of galaxies. An approximate solution of the vacuum 



field equations on the brane, corresponding to weak gravitational fields, is also obtained, and the 



. expressions for the dark radiation and dark mass are derived. The qualitative behavior of the dark 

o . 

, mass is similar to that of the observed virial mass in clusters of galaxies. We compare our model 



with the observational data for galaxy clusters, and we express all the physical parameters of the 



^ ■ model in terms of observable quantities. In particular, we predict that the dark mass must extend 

H : 

far beyond the presently considered virial radius. The behavior of the galaxy cluster velocity dis- 
persion in brane world models is also considered. Therefore the study of the matter distribution 
and velocity dispersion at the extragalactic scales could provide an efficient method for testing the 
multi-dimensional physical models. 
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I. INTRODUCTION 

Several recent astrophysical observations [l| have provided the astonishing result that 
around 95 — 96% of the content of the universe is in the form of dark matter + dark energy, 
with only about 4 — 5% being represented by baryonic matter. More intriguing, around 70% 
of the energy- density may be in the form of what is called the dark energy, which provides 
Qde ~ 0.7, and may be responsible for the acceleration of the distant type la supernovae. 
The best candidate for the dark energy is the cosmological constant A, which is usually 
interpreted physically as a vacuum energy, with energy density p\ and pressure p\ satisfying 
the unusual equation of state oa. = —Va/c^ = A/SnG. Its size is of the order A ^ 3 x 10^^^ 

n 

cm~^ [2]. The analysis of the cosmic microwave radiation background data from the WMAP 
experiment indicated that most of the 22% dark matter content of the Universe must be of 
non-baryonic nature. More specifically, the analysis of the power spectrum indicates that 
a theory of gravity based essentially on the properties of baryonic matter would produce a 
lower third peak [3|. The same results can also be obtained from the angular-size - redshift 
relation for ultra-compact radio sources ^. 

The problem of the dark matter is a long standing problem in modern astrophysics. Two 
important observational issues, the behavior of the galactic rotation curves and the mass 
discrepancy in clusters of galaxies led to the necessity of considering the existence of the 
dark matter at a galactic and extra-galactic scale. 

The rotation curves of spiral galaxies are one of the best evidences showing the prob- 
lems Newtonian mechanics and/or standard general relativity has to face on the galac- 
tic/intergalactic scale. In these galaxies neutral hydrogen clouds are observed at large dis- 
tances from the center, much beyond the extent of the luminous matter. Since the clouds 
move in circular orbits with velocity Vtg{r), the orbits are maintained by the balance be- 
tween the centrifugal acceleration v^^/r and the gravitational attraction force GM{r)/r'^ of 
the total mass M(r) contained within the orbit. This allows the expression of the mass 
profile of the galaxy in the form M(r) = rv^g/G. 

Observations show that the rotational velocities increase near the center of the galaxy 

P 

and then remain nearly constant at a value of Vtgoo ~ 200 km/s [5]. This leads to a mass 
profile M(r) = rv^g^/G. Consequently, the mass within a distance r from the center of the 
galaxy increases linearly with r, even at large distances where very little luminous matter 
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can be detected. 

The second astrophysical evidence for dark matter comes from the study of the clusters 
of galaxies. The total mass of a cluster can be estimated in two ways. Knowing the motions 
of its member galaxies, the virial theorem gives one estimate, My, say. The second is 
obtained by separately estimating the mass of each individual members, and summing these 
masses, to give a total baryonic mass M. Almost without exception it is found that My is 
considerably greater than M, My > M, typical values of My jM being about 20-30 [sf. 

These behaviors at a galactic/extra-galactic scale are usually explained by postulating 
the existence of some dark (invisible) matter, distributed in a spherical halo around the 
galaxies. The dark matter is assumed to be a cold, pressureless medium. There are many 
possible candidates for dark matter, the most popular ones being the weekly interacting 
massive particles (WIMP). Their interaction cross section with normal baryonic matter, 
while extremely small, are expected to be non-zero and we may expect to detect them 
directly. It has also been suggested that the dark matter in the Universe might be composed 
of superheavy particles, with mass > 10^'^ GeV. But observational results show the dark 
matter can be composed of superheavy particles only if these interact weakly with normal 
matter or if their mass is above 10^^ GeV jsl. 

Dark matter may also be in the form of Warm Dark Matter particles, consisting either 
of sterile neutrinos, with masses of several keV, or from early decoupled thermal relics 3]- 
However, constraints based on X-ray fluxes from the Andromeda galaxy have shown that 
dark matter particles cannot be sterile neutrinos, unless they are produced by a nonstandard 
mechanism (resonant oscillations, coupling with the inflaton field), or they have been diluted 
by some large entropy release. The X-rays produced by the decays of these relic sterile 
neutrinos can boost the production of molecular hydrogen, which can speed up the cooling 
of gas and the early star formation (rf. 

Hence, despite more than 20 years of intense experimental and observational effort, up to 
now no non- gravitational evidence for dark matter has ever been found: no direct evidence 
of it and no annihilation radiation from it. Moreover, accelerator and reactor experiments do 
not support the physics (beyond the standard model) on which the dark matter hypothesis 
is based upon. 

Therefore, it seems that the possibility that Einstein's (and the Newtonian) gravity breaks 
down at the scale of galaxies cannot be excluded a 'priori. Several theoretical models, based 
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on a modification of Newton's law or of general relativity, have been proposed to explain 
the behavior of the galactic rotation curves [8|. 

The proposal by Randall and Sundrum that our four- dimensional space-time is a three- 
brane, embedded in a five- dimensional space-time (the bulk), had attracted a considerable 
interest lately. According to the brane-world scenario, the physical fields (electromagnetic, 
Yang-Mills etc.) in our four- dimensional Universe are confined to the three brane. These 
fields are assumed to arise as fluctuations of branes in string theories. Only gravity can freely 
propagate in both the brane and bulk space-times, with the gravitational self-couplings not 
significantly modified. Even if the fifth dimension is uncompactified, standard AD gravity 
is reproduced on the brane. Hence this model allows the presence of large, or even infinite 
non-compact extra dimensions. Our brane is identified to a domain wall in a 5-dimensional 
anti-de Sitter space-time. For a review of dynamics and geometry of brane Universes see 

0. 

Due to the correction terms coming from the extra dimensions, significant deviations from 
the Einstein theory occur in brane world models at very high energies Gravity is largely 
modified at the electro-weak scale 1 TeV. The cosmological and astrophysical implications 



of the brane world theories have been extensively investigated in the physical literature |12 |. 

The vacuum field equations on the brane have been reduced to a system of two ordinary 
differential equations, which describe all the geometric properties of the vacuum as functions 



th e projections of the Weyl curvature of the 
Several classes of vacuum solutions of 



of the dark pressure and dark radiation terms 
bulk, generating non-local brane stresses) in 
the static gravitational field equations in the brane world scenario describing conformally 
symmetric gravitational fields in the vacuum have been obtained in 1^, [isl . As a possible 
physical application of these solutions the behavior of the angular velocity Vtg of the test 
particles in stable circular orbits on the brane has been considered. In the obtained models in 
the limit of large radial distances and for a particular set of values of the integration constants 
the angular velocity tends to a constant value. This behavior is typical for massive particles 
(hydrogen clouds) outside galaxies and is usually explained by postulating the existence 
of the dark matter. Thus, the rotational galactic curves can be naturally explained in brane 
world models without introducing any additional hypothesis. The galaxy is embedded in 
a modified, spherically symmetric geometry, generated by the non-zero contribution of the 
Weyl tensor from the bulk. The extra-terms, which can be described in terms of the dark 
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radiation term U and the dark pressure term P, act as a "matter" distribution outside the 



galaxy 
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17|. 



The existence of the dark radiation term generates an equivalent mass term Mu, which 
is linearly increasing with the distance, and proportional to the baryonic mass of the galaxy 
Mbi Mjj{r) MB{r/rQ). The particles moving in this geometry feel the gravitational 
effects of U, which can be expressed in terms of an equivalent mass. In the framework of 
this model all the relevant physical parameters (metric tensor components, dark radiation 
and dark pressure terms) can be obtained as functions of the tangential velocity, and hence 
they can be determined observationally. Using the smallness of the rotational velocity, 
a perturbative scheme for reconstructing the metric in a galactic halo in the braneworld 

n 

models with induced gravity was developed in [18|. In the leading order of expansion, at 
sufficiently large distances, the obtained result reproduce that obtained in the Randall- 
Sundrum braneworld model 
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17|. 



Similar interpretations of the dark matter as bulk effects have been also considered in 

H. 

The exact galactic metric, the dark radiation and the dark pressure in the fiat rotation 
curves region in the brane world scenario has been obtained in |l^. The deflection of 
photons has also been considered and the bending angle of light is computed. The bending 
angle predicted by the brane world models is much larger than that predicted by standard 
general relativistic and dark matter models. Therefore the study of the light deflection by 
galaxies and the gravitational lensing could discriminate between the different dynamical 
laws proposed to model the motion of particles at the galactic level and the standard dark 
matter models. 

Because of its generality and wide range of applications, the virial theorem plays an im- 
portant role in astrophysics. Assuming steady state, one of the important results which can 
be obtained with the use of the virial theorem is to deduce the mean density of astrophysical 
objects such as galaxies, clusters and super clusters, by observing the velocities of test par- 
ticles rotating around them. Hence the virial theorem can be used to predict the total mass 
of the clusters of galaxies. The virial theorem is also a powerful tool for stability studies. In 
a general relativistic framework several versions of the virial theorem have been proposed 



20l | . Versions of the virial theorem including the effect of a cosmological constat have been 



derived in [2l|] and [221], respectively. 
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It is the purpose of the present paper to consider the virial theorem in the framework of the 
brane world models. By taking into account the effects of the non-compact extra-dimensions 
we derive, by using the coUisionless Boltzmann equation, a generalized virial equality, which 
explicitly takes into account the presence of the bulk effects. In the case of static spherically 
symmetric systems the transmitted projections of the bulk Weyl tensor can be described in 
terms of two quantities, called the dark radiation and the dark pressure, respectively. The 
dark radiation term gives an effective contribution to the gravitational energy, the total virial 
mass being proportional to the mass associated to the effective mass of the dark radiation. 
This term may account for the well-known virial theorem mass discrepancy in actual clusters 
of galaxies. The generalized virial theorem in the brane world models can be an efficient 
tool in observationally testing the brane world models. 

An approximate solution of the vacuum field equations on the brane, corresponding to 
weak gravitational fields, is also obtained, and the expressions for the dark radiation and 
dark mass are derived. The qualitative behavior of the dark mass is similar to that of the 
observed virial mass in clusters of galaxies. 

We consider several astrophysical applications of the general results obtained in the 
present paper. The expressions of the dark radiation and of the dark mass can be ex- 
pressed in terms of observable astrophysical quantities, like the temperature and the gas 
density profile, which can be obtained from the X-ray observations of the gas in the cluster. 
Prom the obtained expression of the dark radius, which describes the relevant length scale 
for brane world effects, we predict that the dark mass must extend far beyond the presently 
considered virial radius. The behavior of the galaxy cluster velocity dispersion in brane 
world models is also considered in detail. Therefore the study of the matter distribution and 
velocity dispersion at the extragalactic scales could provide an efficient method for testing 
the multi-dimensional physical models. 

The present paper is organized as follows. The basic equations for static spherically 
symmetric gravitational fields on the brane are presented in Section II. The generalized virial 
theorem, including the effects of the non-compact extra dimensions is derived in Section III. 
In Section IV we derive the general expressions of the cluster metric, of the dark radiation and 
of the dark mass for the clusters of galaxies. Some astrophysical applications are presented 
in Section V. We discuss and conclude our results in Section VI. 
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II. STATIC SPHERICALLY SYMMETRIC GRAVITATIONAL FIELDS ON THE 
BRANE 



On the 5-dimensional space-time (the bulk), with the negative vacuum energy A5 and 
brane energy-momentum as source of the gravitational field, the Einstein field equations are 
given by 

Gij = klTjj, Tu = -A,gu + 5(F) [-\,grj + T,7"-] , (1) 

with Aft the vacuum energy on the brane and /c| = SttG^. In this space-time a brane is 
a fixed point of the Z2 symmetry. In the following capital Latin indices run in the range 
0, ...,4, while Greek indices take the values 0, ...,3. 

Assuming a metric of the form c/s^ = {ninj+gij)dx^dx'^ , with njdx^ = dx the unit normal 
to the X =constant hypersurfaces and gjj the induced metric on x =constant hyp ersurf aces, 
the effective four-dimensional gravitational equations on the brane take the form 

Gf^i, = —A.g^y + k^T^iy + k^Sp^p — E^u-, (2) 

where S^y is the local quadratic energy-momentum correction 

11 1 

and Efj^iy is the non-local effect from the free bulk gravitational field, the transmitted pro- 
jection of the bulk Weyl tensor Ciajb-, Ejj = CiAjB^^n^ ■, with the property Eu — >■ 
E^uSj6j as X 0. We have also denoted kf = SnG, with G the usual four- dimensional 
gravitational constant. Eq. ([2]) is also known as the effective Einstein equation. It has 
been shown that for a large class of generalized Randall- Sundrum type models, the charac- 
terization of brane-gravity sector by the effective Einstein equation, Codazzi equation and 
the twice-contracted Gauss equation is equivalent with the bulk five-dimensional Einstein 



equation 



23| 



The four-dimensional cosmological constant. A, and the four-dimensional coupling con- 
stant, A;4, are given by A = /c^ (A5 + fc^A^/G) /2 and A;| = k^Xh/6, respectively. The four- 
dimensional gravitational constant G is given by G = k^Xb/ASir 111 ]. 

The Einstein equation in the bulk and the Codazzi equation also imply the conservation of 
the energy-momentum tensor of the matter on the brane, DyT^'^ = 0, where denotes the 
brane covariant derivative. Moreover, from the contracted Bianchi identities on the brane it 



follows that the projected Weyl tensor should obey the constraint D^E^'^ = klD^Su" 



io|. 



In the limit A5 ^ ^ we recover standard general relativity. An alternative possibility in 
recovering the four-dimensional Einstein equations is to take the limit ^5 0, while keeping 



the Newtonian gravitational constant G finite [llj. As one can see from its definition, 
the existence of the four-dimensional gravitational constant relies on the existence of the 
vacuum energy Af,. It is impossible to define an the gravitational constant in an era when 
the distinction between vacuum energy and matter is ambiguous. Moreover, the positivity 
of G also requires the positivity of A;, 

The symmetry properties of E^j^y imply that in general we can decompose it irreducibly 
with respect to a chosen 4-velocity field as 



E 



U ( U^U^ + - V ) + '^QipUu) + P^,u 



(4) 



where k 
term U 



^5/^4, h/j^iy = g^u + u^Ui, projects orthogonal to m'^, the "dark radiation" 
—k~^E^yU^u^ is a scalar, = k~'^h'^Eapu^ a spatial vector and P^^ 



h{p "hi,) ^ — \h^uh"^ Eai3 a spatial, symmetric and trace-free tensor |lO |. 
We assume that the matter on the brane consists of an anisotropic fiuid, characterized by 
an effective energy- density pe// 7^ 0, a radial pressure p^^j and a tangential pressure Pe//> 
respectively. Generally Pgjj 7^ Pe//; but for isotropic systems p^jj^ = pi'f}- Hence T^jy 7^ 0, 
and consequently also S^u 7^ 0. E^i, satisfies the constraint D uE^ y = k^D^^S^iy, where 
is the projection (orthogonal to u'^) of the covariant derivative Uj. In an inertial frame at 



any point on the brane we have = 6q and h^^ =diag(0, 1, 1, 1). In the static spherically 



symmetric case Q/, 



pi/ 

and the constraint for E^iy takes the form 



2M 



Id,u + Iua, + d^p,. + a^p,. 



(5) 



where = u^DyU^ is the 4-acceleration. 

In the static spherically symmetric case we may choose = A{r)r^ and P^v = 
P{r) (r^r,y — where A{r) and P{r) (the "dark pressure" ) are some scalar functions 

of the radial distance r, and is a unit radial vector [25]. 

We choose the static spherically symmetric metric on the brane in the form 



ds' = -e^^-Ut' + e^^'-^dr' + [dO^ + sin^ edct>'') . 



(6) 



Then the gravitational field equations and the effective energy-momentum tensor con- 
servation equation for an anisotropic static spherically symmetric system take the form 
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Q,Q,H,Q,I24| 



+ T + ^ - ^) = 16-Gpi)»+^p.„ (p,„ + 2pit})+^ ([/ - P)-2A. 

(9) 

2C/ + P r(2C/ + P) 2C/ + F ^^^^ 

In the following we shall denote a = IGnG / kfXb = l/^nGXf,. 



III. THE VIRIAL THEOREM IN THE BRANE WORLD MODELS 



We consider an isolated, spherically symmetric cluster, situated in a space with metric 
given by Eq. ([6]). We describe the galaxies, which are treated as identical, coUisionless 
point particles, by a distribution function which obeys the general relativistic Boltzmann 
equation. 

Consider a time-oriented Lorentzian four- dimensional space-time manifold M, with met- 
ric g of signature (— , -|-, -|-, +). The tangent bundle T{M) is a real vector bundle whose 
fibers at a point x G M is given by the tangent space (M). In the space-time M the 
instantaneous state of a particle with mass mo is given by a four-momentum p G Tx (M) at 
an event x G M. The one-particle phase space Pphase is a subset of the tangent bundle given 



by m 



Pphase := {{x,p)\x e M,pe Tx (M) ,p^ = -ml}. (11) 

A state of a multi-particle system is described by a continuous, non-negative distribution 
function / defined on Pphase, and which gives the number dN of the particles of the 

system which cross a certain space-like volume dV at x, and whose 4-momenta p lie within a 
corresponding three-surface element dp in the momentum space. The mean value of / gives 
the average number of occupied particle states {x,p). Macroscopic, observable quantities 



can be defined as moments of / 26|. 



Let {x"}, a = 0, 1, 2, 3 be a local coordinate system in M, defined in some open set 
U C M. The coordinates are chosen so that dt is time-like future directed and da, a = 1,2, 3, 
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are spacelike. Then {d/dx""} is the corresponding natural basis for tangent vectors. We 
express each tangent vector p in f/ in terms of this basis as p = p°'d/dx°' and define a 
system of local coordinates l-^^"^}) A = 0, 7 in Tu{M) as = 2;°+^ = p". This defines 
a natural basis in the tangent space given by [d/dz^] = {a/(9x",9/(9p°} y. 

A vertical vector field over TM is given by tt = p°'d/dp°'. The geodesic flow field a, which 
can be constructed over the tangent bundle, is defined as a = p^d/dx" — p'^p'^V^^^d /dp^ = 
p"Da, where Da = d / dx"' — p'^T'^^d / dp'^ , and are the connection coefficients. Physically, 
a describes the phase flow for a stream of particles whose motion through space-time is 
geodesic [26 1. 

Therefore the transport equation for the propagation of a particle in a curved arbitrary 



Riemannian space-time is given by the Boltzmann equation 

f . d 



26| 



p 



dx'^ 



0. 



(12) 



For many applications it is convenient to introduce an appropriate orthonormal frame or 
tetrad e^"-* (x), a = 0, 1,2,3, which varies smoothly over some coordinates neighborhood U 

•T]"-^ for all X G ?7 



and satisfies the condition e|f^ (x) e|f^ (x 



21 



261]. Any tangent vector p^^ 



at X can be expressed as p^ = p°'e.^a)-> "which defines the tetrad components p"". 

In the case of the spherically symmetric line element given by Eq. ([6]) we introduce the 
following frame of orthonormal vectors 2l|, |26[ |: 



■1 p(2) 



(13) 



Let be the four- velocity of a typical galaxy, satisfying the condition u'^u^ = — 1, with 
tetrad components u^'^^ = u^e^^\ The relativistic Boltzmann equation in tetrad components 



IS 



u 



+ '(6){c)" " 



df 



where the distribution function f = f (x'^, u^'^^^ and l[h){c) ~ ^l^l^(b)^(c) Ricci rotation 

coefficients 



21 



26|. By denoting 

(0) 



Ut,U 



(1) 



Ur, U 



(2) 



(3) 



(15) 



and by assuming that the only coordinate dependence of the distribution function is upon 
the radial coordinate r, Eq. f[T^ becomes 21] 



df (I ^dv ul^u\\ df 1 / 9/ , df 



dUr 



-Ur ue-r. 
r \ dufi 



U4> 



duA 



10 



-e^^'^Ufj, cot I 



ue 



dl_ 



dl_ 

duf) 



0. 



(16) 



The spherical symmetry requires that the coefficient of cot 6 be zero, which imphes that / 
is a function of r, Ur and + only. We multiply now Eq. ( |T6l) by murdu, where m is the 
mass of each galaxy and du = durduedu^/ut is the invariant volume element of the velocity 
space. By integrating over the velocity space and by assuming that / vanishes sufficiently 
rapidly as the velocities tend to ±oo, we obtain 



d_ 
dr 



piu. 



du 

dr ^ 



2(u. 



0, 



(17) 



where, at each point, {uD is the average value of ul etc. and p is the mass density. 
By multiplying Eq. (fT7I) by 47rr^ and integrating over the cluster we obtain 21 1 



- I Aixp \{u^) + luf) + lul 





In terms o 
written as 



r dr + 



1 

2 Jo 



R 



u. 



+ iu 



du 
dr 



dr = 0. 



:i8) 



26| 



the distribution function the energy-momentum tensor of the matter can be 

fmu^Uydu, (19) 



T. 



which gives 



Peff = P (Ut) , Peff = P (ul) , Piff = P {uf) = P 



(20) 



Then, with the use of this form of the energy-momentum tensor, by adding the gravita- 
tional field equations Eqs. ©-Q, we immediately obtain the following relation: 

,/2 „/ ^/^/\ 



u 



u 



2 ^ 4 ^ r 



3aU - A + AnGp 



+ (U,.) + (Uf,) + (u 



+ 



47rG 



(21) 



It is convenient to introduce some approximations at this moment. First of all, we assume 
that u and A are small, so that in Eq. (12T1) the quadratic terms can be neglected. Secondly, 
we assume that the galaxies have velocities much smaller than the velocity of the light, so 
that (m^) , {uf) , (m^) << (m? ) ~ 1- Since for clusters of galaxies the ratio of the matter 
density and of the brane tension is much smaller than 1, p/\b « 1, we can neglect in 
Eq. f l2T|) the quadratic term in the matter density. These conditions certainly apply to test 
particles in stable circular motion around galaxies, and to the galactic clusters. 
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Therefore Eqs. (fTSll and (12T!) become 

I ^ dv 



respectively, where 



-2K + - f ATcr^p^dr = 0, (22) 
2 Jo or 



A- = /" 2»p [(«.?) + U) + ('4)1 r=*, (24) 



is the total kinetic energy of the galaxies. The gravitational potential energy Q of the system 
is defined by 

a=-r^dM(r), (25) 
Jo r 

where M(r) is the mass out to radius r, so that dM{r) = Airpr'^dr. The total mass of the 
system is given by M = /(f dM{r) = Att pr^dr. The main contribution to M is represented 
by the baryonic mass of the intra-cluster gas and of the stars, but other particles, like, for 
example, massive neutrinos, may also give a significant contribution to M. 
Multiplying Eq. (!23|) by and integrating from to r we obtain 

1 f)i, 1 rr 

GM{r) = -r^— + -Ar^ - 3a J U{rydr. (26) 
In the following we denote 

GMu (r) = 3a T U{rydr. (27) 
Jo 

This quantity may be called as the dark mass. By multiplying Eq. (126!) with dM{r), inte- 
grating from to R, by introducing the moment of inertia of the system as / = r'^dM[r), 
and by denoting 

Jo r 

we obtain 



I rR Qy 1 

-n= - Airr^p—dr + -A/ - Qu- (29) 
2J0 dr 3 ^ ^ 

Finally, with the use of the Eq. fl22|) . we obtain the generalization of the virial theorem 
in the brane world models in the form 

2K + Q + -AI -Qu = 0. (30) 

3 
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The generalized virial theorem, given by Eq. (130|) . can be written in a simpler form if we 
introduce the radii Ry, Ri and Ru defined by 

Ry = mV r ^^^dM{r), (31) 
JO r 

1/2 



Ri 







R \ 

r^dM{r) /M{r 



(32) 



R^^M^jTMEmm^ (33) 







so that 



(34) 

/ = MR], (35) 
^ f^. (36, 
flf/ may be called the dark radius of the cluster of galaxies. The virial mass My is defined 



as 



2l| 



2A- ^ (37) 

If this expression is substituted into the virial theorem, given by Eq. (l30l) . we obtain 

My _ ( M^Ry _ J^V'' 

where p = SM/ATrRyRj. 

For U = 0, Mu = and we reobtain the virial theorem in the presence of a cosmological 



constant 



2l|. 



If My/M > 3, a condition which is true for most of the clusters, the term unity can be 
neglected in the bracket in Eq. ( l38l) with little loss of accuracy. Moreover, the contribution 
of the cosmological constant to the mass energy of the galaxy can also be ignored as being 
several orders of magnitude smaller than the observed masses. Therefore the virial mass in 
the brane world models is given by 



My^Mu^^. (39) 

From observational point of view the virial mass My is determined from the study of 
the velocity dispersion of the stars and of the galaxies in the clusters. According to our 
interpretation, most of the mass in a cluster with mass M^ot should be in the form of the dark 

13 



mass M(7, so that Mjj Mtot- A possibility of detecting the presence of the dark mass and 
of the astrophysical effects of the extra dimensions is through gravitational lensing, which 
can provide direct evidence of the mass distribution and of the gravitational effects even at 
distances extending far beyond of the virial radius of the cluster. 



IV. DARK MASS, DARK RADIATION AND METRIC FOR CLUSTERS OF 
GALAXIES ON THE BRANE 

The total mass of galaxy clusters ranges from 10^^ Mq for groups up to a few 10^^ Mq 



for very rich systems. The c^ 
peaked main component 28, 



uster morphology is usually dominated by a regular centrally 



29|. 



As clusters are "dark matter" dominated objects, their formation and evolution is driven 
by gravity. The mass function of the clusters is determined by the initial conditions of the 
mass distribution set in the early universe. The evolution of the large scale matter distri- 
bution on scales comparable to the size of the clusters is linear. The overall process of the 
gravitational growth of the density fluctuations and the development of gravitational insta- 
bilities leading to clustei^^formation has been extensively studied by using both analytical 
and numerical methods [30!]. 

An other important component is the intra-cluster gas, which is assumed to be isothermal, 
and in hydrostatic equilibrium. Most of the baryonic mass in the cluster is contained in 
the gas. The properties of the gas are important observational quantities, since the total 
gravitational mass of the cluster Mtot is determined by using the radial gas distribution 



281]. By assuming that the density and pressure of the gas are pg and Pg, respectively, and 
by neglecting the quadratic terms in the field equations, the basic equations describing the 
structure of a galactic cluster in the brane world model are given by 

-e-^l^-jj+^ = 8nGpg + 3aU + A, (40) 



-A 



- + \]-\ = 8^Gpg + a{U + 2P)- A, (41) 



e-'l + Y + - ^) = ^^^^^ + « (f/ - P) - A, (42) 

U' + 2P' 6P ^-^ Pa+Pa ' 
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Eq. (HOj) can immediately be integrated to give 

,-A_i C 2GM, GMujr) A^,^ 
r r r 3 ' 

where C is an arbitrary constant of integration, Mu is defined according to Eq. (127|) and Mp 

is the total mass of the gas. 

By substituting z/' given by Eq. ( H3|) into Eq. ( HTi) and with the use of Eq. ( I44l) we obtain 

the following system of differential equations satisfied by the gas density, pressure and the 

dark radiation term U, the dark pressure P and the dark mass M^/, respectively, describing 

the gravitational field of a cluster of galaxies in the brane world model: 

= ^^Pa^'- (45) 

-^ = ^rU. (46) 

dU_ _ {2U + P){C + 2GMg + GMu + [a{U + 2P) + 8npg] r^} - fAr^ _ 

'dr ~ (l C 2GMg GMu(r) A ^2 

\ r r r 3 

dr r ^ ' 2U + P' 



2^ - — + 2 (4vrG) ^fy^p',. (47) 



Since the most important contribution to the structure and dynamics of the clusters of 
galaxies comes from the "dark matter", which we interpret as a multi-dimensional effect, 
described by the contribution of the Weyl tensor from the bulk, in discussing the properties 
of the clusters we can neglect, in the first approximation, the effect of the ordinary matter, 
by taking Pg ~ and ~ 0, respectively. 

Therefore the metric components z/(r) and A(r) inside the spherically symmetric static 
cluster are related to the dark radiation U and to the dark pressure P via the Einstein grav- 
itational field equations and the effective energy-momentum tensor conservation equation, 
which in the vacuum take the form [13I, [l^ 

^ = '-^r^U. (48) 
dr G 

dU _ {2U + P){G + GMu + [a{U + 2P)] r^} - fAr^ _ dP_6P 

dr ~ r^(i-C_ GMu(r) _ K^2\ dr r ' ^ ' 

\ r r 3 / 

The system of equations ( 148|) and ( H9|) can be transformed to an autonomous system of 
differential equations by means of the transformations 

C^GMc/^Aa n = Sar^U + 3r^A, p = Sar"^ P - 3r^ A, e = \nr. (50) 
r r 3 
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We shall call and p the "reduced" dark radiation and pressure, respectively. 
With the use of the new variables given by Eqs. fl50l) . the system of Eqs. psj) and fH9|) 
become 

d9=^-'^ (51) 

df, (2/i + p)[g + |(/i + 2p)] dp 

2— + 2^-2p, (52) 



de 1-q de 

In order to close the system of equations Eqs. (15T|) and (!52|) . we need to specify the 
"equation of state" of the dark pressure. As a possible equation of state relating p and /x we 
will assume a general linear relation of the form 

p(/x) = (r-2)/x + S, (53) 

with r and B arbitrary constants. Therefore Eq. (|52l) takes the form 

f^l(2T-3),] = JIt±m±(^:^:IMl±im^2(3-T),-2B, (64) 

where we have neglected the possible effect of the cosmological constant on the structure of 
the cluster. 

For clusters of galaxies with virial masses of the order of My ^ lO^'^M© and virial radii of 
the order of Ry ~ 2 Mpc, the quantity GM/R is of the order of 2 x 10^^. Since observations 
show that inside the cluster the mass is a linearly increasing function of the radius r, the value 
of this ratio is roughly the same at all points in the cluster. Therefore from its definition it 
follows that generally g << 1, and 1 — g ^ 1. Moreover, the quantities and qdq/de are also 
very small as compared to q. Eq. (15T1) gives fi = q + dq/de, dfi/de = dq/de + d'^q/de^. Hence, 
by neglecting the second order terms, we obtain for q the following differential equation: 

^ + m^ + ng = 6, (55) 



where we denoted 



and 



^ B 2r 5 + 3 -18 

m = l-\ \ , (56) 

3 3 (2r - 3) ' ^ ^ 

B 2(i? + 3)r + 3(i?-6) 
''-3"+ 3 (2r - 3) ' ^^^^ 

_ 2B{B + 3) 

^ - 3(3-2r)' ^^^^ 
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respectively. The general solution of Eq. (l55i) is given by 

g(e) = go + Cie'^' + C2e'^', (59) 
where Ci and C2 are arbitrary constants of integration, and we denoted 

go = -, (60) 

n 

and 

^1,2 = ^ , (61) 

respectively. The reduced dark radiation term is given by 

/i (6) = go + Ci (1 + h) e''' + C2 (1 + I2) e'''. (62) 



For a positive m from Eq. ( 1611) it follows that both li and I2 are negative numbers, /i < 
and /2 < 0, respectively. In the original radial variable r we obtain for the dark radiation 
and the mass distribution inside the cluster the expressions 



3aU{r) = ^ + Ci (1 + k) r'^-^ + C2 (1 + k) r''''\ (63) 
and 

GMxjir) = gor I 1 + — r'^ + — r'^ - — ) , (64) 

V go go q^r ) 

respectively, where we have neglected the possible effect of the cosmological constant. 
The metric coefficient v can be calculated from the equation 

^'(^) = g(^) + ^^/^(^^) + ^, (65) 



givmg 



■ 3 + (2r-3)(l + /i) ^, 3 + (2r-3)(l + /2) ^,, 
^ 3/1 ^ 3/2 



(66) 



where C3 is an arbitrary integration constant. In the limit of large distances e^'-''^ behaves 
like e'^*^''^ ~ (7gj"2{(jor+B)/3^ Inside the cluster we can approximate e~'^ ~ 1 — C/r — GMu{r)/r. 

The first term in the mass profile of the dark mass given by Eq. fl64l) is linearly increasing 
with r, thus having a similar behavior to the dark matter in clusters of galaxies. In order to 
obtain a model consistent with observations, the second and third terms in Eq. (!64l) must 
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be decreasing with increasing r. Therefore the constants h and I2 must also satisfy the 
conditions /i + 1 < and Z2 + 1 < 0, respectively. By assuming that m and n are positive 
numbers, the conditions li < and ^2 < are realized automatically. 

In the particular case i? = 0, we have go = 0, and the solution is defined only for values of 
r so that r > 21/8. This shows that a radiation like equation of state of the form P = U/3 
is not allowed in the present model. Moreover, the obtained mass profile is not consistent 
with the observations, and for such an equation of state the dark radiation cannot play the 
role of the dark matter. For the case T = 4, corresponding to an equation of state P = 2U, 
the dark mass is given by GMu{r) = Ci/r^'^'^ + C2/r^'^'^ — C. 

Since according to our physical interpretation Mjj is an effective, geometry induced mass, 
it must satisfy the condition > for all r. Therefore the solution obtained in the present 
Section is valid only for values of the coordinate radius r so that gor + Cir'i+^ + C2r'2+^ — C* > 
0, where we assume that all integration constants are non-zero. In the limit of small r, taking 
into account that /i + 1 < and ^2 + 1 < 0, and by assuming that the constant C is small 
and can be neglected near the origin of the cluster, we obtain Cir'-^ + C2r'^ ^ 0. If both 
Ci and C2 are strictly positive, then the dark mass diverges at the center of the cluster, 
lim^^o Mu{r) = 00. On the other hand, if Ci and C2 have different signs, there is a single 
point ro, given by ro ~ (C2/Ci)^^'''^~'^\ so that Mu{ro) ^ and Mu{ro) < 0, for r < tq. In 
this case our solution has physically acceptable properties only in the region r > tq. 

The solution of the brane world model Einstein field equations, given by Eqs. ( |64l) . (|63l) 
and [66]) depends on five, equation of state related or integration constants, go, C and Ci, 
i = 1,2,3, respectively. These constants must be determined by using some appropriate 
boundary conditions, or by using some other observable quantities. Without any loss in 
generality we can take the integration constant C3 = 1, since the metric tensor component 
exp (z/) can be arbitrarily re-parameterized by a transformation of the time coordinate. 

An important observational quantity is the radial velocity dispersion a^, which is related 



to the total mass in the cluster by the relation GMy = cr^Rv- 271]. By assuming that 
Ru ~ Rv gives, with the use of the generalized virial theorem. My (Ry) ~ Mu (Ry), 
respectively. For r = Ry Eq. flMl) becomes Mu (Ry) ~ %Ry, thus giving 

go ~ c^l- (67) 
Observationally, the mass profiles of the clusters of galaxies are obtained from the equa- 
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tion 



GMtot{r) = -aUr)r 



d\n cr;{r) 



+ 



dlnnir] 



(68) 



dlnr dlnr 

where n(r) is the spatial galaxy number density and f3cr is the velocity anisotropy parameter 
271]. This equation has a very similar mathematical structure with Eq. (IMl) . with the total 
observed mass of the cluster interpreted as the dark mass. Therefore the other constants in 
the model can be, at least in principle, obtained by fitting the dark radiation U{r) and dark 
mass Mij{r) with the observed density and mass profiles of the dark matter in clusters of 
galaxies. 



V. ASTROPHYSICAL APPLICATIONS 



Astrophysical observations together with cosmological simulations have shown that the 
virialized part of the cluster corresponds roughly to a fixed density contrast 6 ~ 200 as 
compared to the critical density of the universe, Pc{z), at the considered redshift, so that 
Pv = SMy/AttRy = 6pc (z), where py is the virial density. My and Ry are the virial mass 
and radius, Pc{z) = h'^{z)3HQ/8TcG, and h{z) is the Hubble parameter normalized to its 
local value: h'^{z) = Qm (1 + z)^ + where Q^n is the mass density parameter and is 
the dark energy density parameter ^]. 

Once the integrated mass as a function of radius is determined for galaxy clusters, a 
physically meaningful fiducial radius for the mass measurement has to be defined. The radii 
commonly used are either r2oo or '"soo- These radii are the radii within the mean gravitational 
mass density of the matter (ptot) = 200pc or SOOpc- A pragmatic approach to the virial mass 
is to use r2oo as the outer boundary [28|. 

The numerical values of the radius r2oo are in the range r2oo = 0.85 Mpc (for the cluster 
NGC 4636) and raoo = 4.49 Mpc (for the cluster A2163). A typical value for r2oo is 2 Mpc. 
The masses corresponding to r2oo and rsoo are denoted by M200 and M500, respectively. 
Usually it is assumed that My = M200 and Ry = T200 



A. Dark radiation, dark mass and dark radius from galactic cluster observations 

In clusters of galaxies most of the baryonic mass is in the form of the intra-cluster gas. 
The gas mass density Pg distribution can be fitted with the observational data by using the 
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following expression for the radial baryonic mass (gas) distribution 28 1 

,2\ -3/3/2 
I2 



p,(r)=po 1 + ^1 , (69) 



where Tc is the core radius, and po and [3 are (cluster-dependent) constants. 

A static spherical system in collisionless equilibrium can be described by the Jeans equa- 
tion, which is given by 

d 



dr 



+ ^ {-'r - 4.) = (0 % (70) 

where $ (r) is the gravitational potential, ar and ue^^ are the mass-weighted velocity dis- 
persions in the radial and tangential directions, respectively. We assume that the gas is 
distributed isotropically inside the clusters and therefore we take cr^ = ae^^. The pressure 
profile Pg is related to the velocity dispersion by Pg = PgCr"^. 

Therefore the equation describing the isotropic equilibrium of the gas is 



dPg (r) _ ^ d<!> (r) _ GMtot (r 
dr dr 



P. (r) = -^^^^P. (r) , (71) 



where Mtotij) is the total mass inside radius r. In obtaining Eq. (ITTi) we have assumed that 
the gravitational field is weak, and that the gravitational potential satisfies the usual Poisson 
equation A$ ^ 4:7rptot, with the total energy density ptot ~ Pg+Pm + 3«f//87rG also including 
the energy density Pm of other forms of matter, different from gas, like, for example, luminous 
matter, massive neutrinos etc., and the bulk (five-dimensional) effects. In the following we 
will systematically neglect the quadratic contribution to the energy-momentum tensor. The 
energy density of the gas is much larger than the pressure, Pg >> Pg. 

The observed X-ray emission from the hot, ionized intra-cluster gas is usually interpreted 
by assuming that the gas is in isothermal equilibrium. Therefore, we may assume that the 
pressure Pg of the gas satisfies the equation of state Pg (r) = (kp Tg/ urrip) pg (r), where fc^ is 
Boltzmann's constant, Tg is the gas temperature, p ~ 0.61 [28| is the mean atomic weight 
of the particles in the cluster gas, and rup is the proton mass. Then Eq. (ITTil gives 



By using the density profile o: 
inside the cluster the relation 



the gas given by Eq. (!69l) we obtain for the mass profile 



f . ^kePTg r^ 
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On the other hand, from the Einstein field equation Eq. ([7j) it follows that the total mass 
inside the radius r satisfies the following mass continuity equation, 

dMtot (r) 



dr 



Anr pg (r) + 3ar U (r) 



(74) 



where we have neglected again the quadratic contributions in the gas density and pressure 
to the total mass. Since the gas density and the total mass profile inside the cluster are 
given by Eqs. (1691) and ( 1731) . respectively, we can obtain immediately the expression of the 
dark radiation term inside the cluster as 



3aGU{r) 



1 



(75) 



In the limit r >> Tc we obtain for the dark radiation the simple relation 



3aGU{r) 



^-MIl-ATiGp^rfr^-^^ 
prUp 



1 



(76) 



The dark mass can be obtained generally as 



GMuir) = 3a f r'^U{r)dr 
Jo 



AnGpo 



r'^dr 



3/3/2 



prup 1 + r2/r2 ' " Jo (l + r2 /r2) 

In the limit r » Tc the dark mass Mu is given as a function of r by the expression 



(77) 



GMu (r) 



r. 



(78) 



pnip 3 (1 — /5) 

Let's assume first that we can neglect the contribution of the gas to the dark radiation 
and dark mass, respectively. Then we obtain 

[3kBf3Tg\ _2 



and 



3aGU{r) 



GMu (r) 



prUp ) 

'"iks^Tg^ 
, /imp ^ 



(79) 



(80) 



respectively. 

There is a simple way to estimate an upper bound for the cutoff of the dark mass. The 
idea is to consider the point at which the decaying density profile of the dark radiation 
associated to the galaxy cluster becomes smaller than the average energy density of the 
Universe. Let the value of the coordinate radius at the point where the two densities are 

Puniv, where puniv is the mean energy 



equal to be J^^^. Then at this point 3aGlJ 



(cr) 
V 
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density of the universe. By assuming puniv = Pc = SH'^/SnG = 4.6975 x 10 ^^/^-so g/cm ^, 
where H = 50/z5okm/Mpc/s [28j, we obtain 



The total dark mass corresponding to this value is 

M^^ = Mu = 4.83 X lO^V'/^ h^o'Me. (82) 

This value of the mass is consistent with the observations of the mass distribution in the 
clusters of galaxies. However, according to the brane world scenario, we predict that the 
dark mass and its effects extends beyond the virial radius of the clusters, which is of the 
order of only a few Mpc. 

Let's now briefly consider the effect of the gas on the dark mass distribution. If /5 < 2/3, 
the dark mass has a maximum value, which occurs at a radius -Ri, which can be obtained 
from the condition dMu{r) / dr = 0, and is given by 



The maximum value of the dark mass can be found as 

GMjr (flO = gf^ (4.Gpy/r""'> . (84) 



3 (1 — /3) Y pTUp 

The central density of the gas po can be estimated from the total gas mass Mg as 

3{l-(3)Mg{Rv) 



Hence we obtain 

Ri 

and 



kepTg 



(1 - p) GMg (Rv) prup 



1/(2-3/3) 



^3(l-/3)/(2-3/3)^ (86) 



3(1-/3) 

CM- {R,) = [GMg (i?,)]-/(-^^) Rf-m2~s,)_ (87) 

[3 (1 — pjj \ fJ'^p J 



Usually the mass of the gas represent 5 — 10% of the total mass in the cluster 27|, |28 |. 
For /3 = 1/2, kBTg = 5 keV, and assuming for the mass of the gas a value of the order 
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Mg = 7 X lO^Mo, we obtain Ri ^ 50.92 Mpc and M^^^ (i?i) = 6.6 x IO^^Mq, respectively. 
For /3 > 2/3 the mass is a linearly increasing function of the radius. 

By comparing Eqs. (15T]) and (15^ we predict that, due to the presence of the gas in the 
cluster, the dark mass, playing the role of the "dark matter", should cut-off earlier, and 
before it merges into the background. 

Astrophysical observations are extended and interpreted in terms of an outer density of 
the order of Spuniv, with a fixed density contrast 6 = 200. In order to compare the predictions 
of the brane world models with the observations it is necessary to estimate the dark radius 
for an outer density of the dark radiation of the order of 6pcr, with the density contrast 
6 = 200. Hence, the dark radius Rfj^ corresponding to the density contrast 6 = 200 is given 
ClS db solution of the non-linear algebraic equation 

^ - 4-Gp„ (^) (ijr) ' = i Gp. (ijr) ' . (88) 

By evaluating the gas density at Rfj^ we obtain 

/ \ 3/3 

^ Po j =Sg,Pc, (89) 

where Sgd gives the density of the amount of the total mass in form of the gas at a distance 
Rij^ from the center of the cluster with respect to the critical density. Then we obtain for 
the dark radius of the cluster the expression 



^200 ^ ^kB/3T, <^ 91.33 rJ kBTX^\^^ 

For ksTg = 5 keV, (3 = 2/3 and 6gd = 20 we obtain i^^oo ^ 3.51 Mpc. This shows that 
is of the same order of magnitude as the virial radius Ry = T2005 and the presence of 
the gas considerably decreases the magnitude of the dark radius. Therefore it follows from 
the virial theorem that the virial mass of the galaxy clusters is mainly determined by the 
dark mass associated to the transmitted projection of the bulk Weyl tensor. 



B. Radial velocity dispersion in galactic clusters 

Th e virial mass can also be expressed in terms of the characteristic velocity dispersion o"i 

Mv = ^alRv. (91) 



as 



27| 
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By assuming that the velocity distribution in the cluster is isotropic, we have {u^) = 
{uf.) + (uq) + (u^^ = 3 (u^) = 3cr^, with cr^ the radial velocity dispersion, ai and cr^ are 
related by Saf = a^. 

In order to derive the radial velocity dispersion relation for clusters of galaxies in brane 
world models we start from Eq. fll7p . Taking into account that the velocity distribution is 
isotropic, we obtain 

Since inside the cluster ~ 1, by neglecting the cosmological constant the Einstein 
field equation Eq. fl2T|) becomes 

v" v' 

_ + _ = 3af/ + 47rG'p, (93) 
2 r 



and can be written in the equivalent form 

1 d 
2r^ dr 



At ^ f r V) = 3af/ + 47rG'p, (94) 
n.r \ ' 



giving upon integration 

r^z/' = 2GMu{r) + 2GM(r) + 2C^, (95) 

where C4 is an arbitrary constant of integration. Since from Eq. fl92|) we have u' = 

— {2/p)d{pal) /dr, it follows that the radial velocity dispersion of the galactic clusters on 

the brane satisfies the differential equation 

d / 2\ GMuir) , , GM(r) , ^ C4 , , 
- (por^) = ^p(r) ^p(r) - 4p(r), 96 

with the general solution given by 



GMu{r) , , GM[r) . , a 



rfr + ^, (97) 



P 

where C5 is an integration constant. 

As an example of the application of Eq. ( 1971) we consider the case in which the density 
p of the normal matter inside the cluster has a power law distribution, so that 

p(r) = p^r-\ (98) 

with Po and 7 7^ 1, 3 positive constants. The corresponding normal matter mass profile is 
Mir) = Airpor^^'' / (3 — 7). The dark mass GMu = q^r is linearly proportional to r, as has 
been shown in the previous Section. Therefore after some simple calculations we obtain 

a^ r = — + 7 TTT^ ^ + ^— - + — r^, 7 ^ 1, 3. (99) 

7 (7 - 1) (3 - 7) 7 + 1 r Po 
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In the case 7 = 1 we obtain 



(^rir) = go + 7T^ - 27iGpor\nr + ^r,7 = 1. (100) 

For 7 = 3 we find 

= f - .Gpo + + + = 3. (101) 

Usually the observed data for the velocity dispersion in clusters of galaxies are analyzed by 
assuming the simple form crf.{r) = B/{r + b) for the radial velocity dispersion, with B and b 
constants. As for the density of the galaxies in the clusters the relation p (r) = A/r (r + a)^, 
with A and a constants, is used. The data are then fitted with these functions by using a 
non-linear fitting procedure For r << a, p{r) ^ A/r, while for r >> a, p{r) behaves 
like p(r) ^ A/r^. Therefore the comparison of the observed velocity dispersion profiles of 
the galaxy clusters and the velocity dispersion profiles predicted by the brane world scenario 
may give a powerful method to discriminate between the different theoretical models. 



VI. DISCUSSIONS AND FINAL REMARKS 



The galactic rotation curves and the mass distribution in clusters of galaxies continue to 
pose a challenge to present day physics. One would like to have a better understanding of 
some of the intriguing phenomena associated with them, like their universality and the very 
good correlation between the amount of dark matter and the luminous matter in the galaxy. 
To explain these observations, the most commonly considered models are based on particle 
jhysics in the framework of Newtonian gravity, or of some extensions of general relativity 

In the present paper we have considered, and further developed, an alternative view to the 

nnnn 

dark matter problem p^, 115|, 116|, , namely, the possibility that the galactic rotation curves 
and the mass discrepancy in clusters of galaxies can naturally be explained in models in 
which our Universe is a domain wall (a brane) in a multi-dimensional space-time (the bulk). 
The extra-terms in the gravitational field equations on the brane induce a supplementary 
gravitational interaction, which can account for the observed behavior of the galactic rotation 
curves. By using the simple observational fact of the constancy of the galactic rotation 
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curves, the galactic metric and the correspon ding Weyl stresses (dark radiation and dark 
pressure) can be completely reconstructed |l6l. llTI]. 

In order to estimate the effect of the bulk effects on the extra-galactic scale, at the level of 
the clusters of galaxies, we have generalized the virial theorem to include the contribution of 
the non-compact extra dimensions. To derive the generalized virial theorem we have used a 
method based on the collisionless Boltzmann equation. The generalized virial theorem allows 
to remove the virial mass discrepancy, by showing that the virial mass My is proportional, 
as one can see from Eq. (!39|) . to the dark mass Mu associated to the dark radiation term 
t/, which is the scalar component of the projected Weyl tensor E^^. The dark mass can be 
directly related to the observed virial mass, and, for typical clusters of galaxies, it is of the 
same order of magnitude as the virial mass. This shows that at the extra-galactic scale the 
dark mass plays the role of what is conventionally called the dark matter. 

In the present model the cut-off of the dark mass (which represents the dark matter in 
the conventional interpretation) is not due to the merging of the dark radiation with the 
cosmological background, but it is due to the presence of the intra-cluster gas, which is the 
only form of conventional matter in our theory. 

In order to derive our main results we have used a set of approximations and assumptions, 
which can be classified as a) brane world model related b) small velocity limit approximation 
and c) weak gravitational field approximation, respectively. In our treatment of the brane 
world models we have systematically ignored the second order terms in the matter density. 
This approximation is justified since the ratio p/A;,, where A;, is the brane tension, is very 
small in the case of the galaxy clusters. By taking for the mean density of the baryonic matter 
in the cluster a value of p ~ 10^^° — 10^^'^ g/cm^ 28|, l29l. l30l| . and for the brane tension a 
value of the order of \ = (100 GeV)^ = 1.37 x 10^^ g/cm^ ^ , it follows that the ratio p/Xb 
is indeed negligible small at the level of the clusters of galaxies, p/\b « 1, and hence the 
second order effects in the matter density can be safely neglected. Since the dispersion of 



the velocity of the galaxies in the clusters is of the order 600 — 1000 km/s 28|, [29!, l30|], giving 
for the square of the ratio of the velocity and of the speed of light a value of the order of 
{v/cy ^ 4 X 10~6 - 1.11 X IQ-^ « 1, we can neglect the relativistic effects in the Boltzmann 
equation, and use the small velocity limit of it. The intensity of the gravitational effects can 
be estimated from the ratio GM/R, which for typical clusters is of the order of 10~^ << 1. 
Therefore inside galactic clusters the gravitational field is weak, and the linear approximation 
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of the gravitational equations provides an accurate description of the dynamics of the massive 
test bodies. We have also systematically neglected the possible effects of the cosmological 
constant. As for the physical structure of the clusters, we have assumed that they are 
dark radiation dominated spherically symmetric static astrophysical systems, with a normal 
baryonic matter content representing only 5 — 10% of the total equivalent mass energy 
density. Both these assumptions are very well supported by observations and cosmological 
simulations Q, Q- 

In the present model all the relevant physical quantities, including the dark mass and 
the dark radius Ru, which describe the non-local effects due to the gravitational field of the 
bulk, are expressed in terms of observable parameters (the temperature and the virial radius 
of the cluster, and the observed velocity dispersion). Hence it follows that the contributions 
of the bulk Weyl tensor in the case of the clusters of galaxies on the brane - the dark mass 
Mjj, the dark radius Ru, and the gravitational energy Qu - can be determined from the 
study of the X-ray properties of the gas within the virial radius. 

Generally, the virial mass My is obtained from the observational study of the velocity 
dispersions of the stars in the cluster. Thus, it cannot give a reliable estimation of the 
numerical value of the total mass Mg + Mjj in the cluster. A very useful method for the 
study of the total mass distribution in clusters is the gravitational lensing of the light, which 
may provide direct evidence for the total mass distribution, and for the gravitational effects 
at large distances from the cluster. 

Therefore, this opens the possibility of testing the brane world models by using astro- 
nomical and astrophysical observations at the extra-galactic scale. In this paper we have 
provided some basic theoretical tools necessary for the in depth comparison of the predictions 
of the brane world model and of the observational results. 
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